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DERIVATION  OF  THE  VARIABLE  COEFFICIENT  KORTEWEG-de  VRIES 
EQUATION  USING  MULTITIMING  TECHNIQUES 


1.  Introduction 

In  the  theory  of  shallow  water  waves  the  Korteweg-de  Vries  (KdV)  equation 

V,  +  j  £  vv  +  -g-  Vxxx  ”  0 

plays  a  fundamental  role.  This  equation  was  first  derived  by  Korteweg  and  de  Vries  (1895)  and 
describes  the  motion  due  to  small  amplitude  waves  in  shallow  water  of  uniform  undisturbed 
depth.  By  shallow  water  we  mean  that  the  depth  of  the  water  is  small  compared  to  some 
characteristic  length  of  the  waves.  The  nondimensional  parameter  £  is  the  Ursell  number  and 
is  defined  by 


where  a  is  the  maximum  amplitude  of  the  wave,  X  is  a  characteristic  length  scale  of  the  wave 
and  h0  is  the  undisturbed  depth  of  the  water.  The  importance  of  the  KdV  equation  is  that  it 
can  describe  slightly  dispersive  linear  wavetrains  (I  «  1),  nonlinear  long  waves  (I  »  1)  or 
cnoidal  and  solitary  waves  I  —  0(1).  Although  the  KdV  equation  cannot  describe  all  the 
phenomena  associated  with  shallow  water  waves,  it  nevertheless  is  an  extremely  useful  tool  in 
the  study  of  such  waves. 

It  is  frequently  of  interest  to  determine  the  behavior  of  these  waves  as  they  move  into 
regions  of  shallower  and  shallower  water.  The  first  study  of  this  prol  'em  is  due  to  Green  (see 
Lamb  (1932),  Sect.  185).  Green  shows  that  when  infinitely  long  linear  waves  travel  over  a  gra¬ 
dually  varying  bottom  the  amplitude  varies  as  )TI/4  where  h  =  h(x)  is  the  undisturbed  depth 
of  the  water.  This  amplification  (or  attenuation)  law  is  often  referred  to  as  Green’s  law.  Other 
authors  (see,  e.g.,  Grimshaw  (1970)  and  Johnson  (1973))  have  studied  the  behavior  of  solitary 
waves  as  they  propagate  into  water  of  decreasing  depth.  There  it  was  shown  that  the  wave 
amplifies  like  h~'.  Svendsen  and  Hansen  (1978)  have  studied  the  behavior  of  cnoidal  waves  as 
they  approach  a  beach  and  determined  the  Green's  law  appropriate  for  these  waves.  In  each  of 
these  cases  the  wave  was  regarded  as  a  slowly  varying  solitary,  cnoidal  or  linear  wave.  In  all  of 
these  studies  but  one,  viz.  Johnson  (1973),*  no  attempt  was  made  to  derive  an  equivalent  KdV 
equation  for  problems  involving  gradual  depth  changes.  Johnson  showed  that  the  correct  equa¬ 
tion  governing  the  surface  elevation  tj  is 

.3,-1  .  i  h*  .  »  <#>  1  „ 

+  t  *"  +  7  m  *  ”  ■ 0 

where  x  is  the  nondimensional  distance  measured  parallel  to  direction  of  propagation,  r  is  a 
wave  coordinate 
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•The  above  equation  for  r)  has  also  be  |iven  without  proof  by  Ostrovskiy  and  Pelinovskiy  (1970). 
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where  t  is  the  nondimensional  time  and  h  =  h(f>x)  is  the  undisturbed  water  depth  seated  with 
some  typical  depth  h0.  Here  «,  /a2,  ft  are  small  parameters  defined  later  which  measure  the 
wave  amplitude,  the  shallowness  of  the  water  and  the  gradualness  of  the  bottom  variation.  It  is 
clear  that  this  equation  is  just  the  flat  bottom  KdV  equation  with  slowly  varying  coefficients. 
Because  of  the  importance  of  this  equation  to  the  study  of  waves  in  water  of  varying  depth  it  is 
the  intent  of  this  paper  to  give  a  second  derivation  of  this  equation  using  the  standard  tech* 
niques  of  singular  perturbation  theory.  The  derivation  of  Johnson  does  not  state  the  limits  on 
the  validity  of  the  above  equation  nor  is  it  clear  for  what  relative  values  of  the  small  parameters 
€,  fi2,  8  the  above  equation  holds.*  Furthermore,  the  derivation  presented  by  Johnson  appears 
to  preclude  the  possibility  of  reflected  waves  at  the  second  order;  see,  for  example,  his  expan¬ 
sion  of  Tj.  Both  higher  order  reflected  and  transmitted  waves  are  known  to  exist.  The  reflected 
waves  have  been  discussed  by  Peregrine  (1967)  and  Madsen  and  Mei  (1969).  In  the  case  of 
long  linear  waves  Witting  (private  communication)  has  shown  that  reflected  and  transmitted 
waves  are  necessary  for  the  satisfaction  of  certain  conservation  laws. 

In  the  following  report  the  variable  coefficient  KdV  equation  is  derived  by  a  straightfor¬ 
ward  application  of  the  method  of  multiple  scales  or  multitiming.  This  technique  has  been  dis¬ 
cussed  by  many  authors;  excellent  references  are  Cole  (1968),  Van  Dyke  (1975),  Nayfeh 
(1974)  and  Leibovich  and  Seebass  (1974).  The  analysis  is  conceptually  straight-forward  and 
because  of  the  widespread  use  of  the  multitiming  technique  it  is  felt  that  the  derivation 
presented  here  is  more  accessible  to  the  reader.  Furthermore,  the  limits  on  the  theory  come 
immediately  out  of  the  analysis.  In  addition,  the  assumptions  made  allow  for  the  possibility  of 
reflected  waves  in  the  higher  order  terms. 

2.  Mathematical  Formulation 

We  consider  two  dimensional  motions  of  an  inviscid  fluid.  The  coordinate  system  is  that 
sketched  in  Figure  1.  The  distance  along  the  undisturbed  free  surface  is  x  and  the  distance 
measured  normal  to  it  is  y.  The  acceleration  of  gravity  is  in  the  negative  y  direction  and  the 

bottom  is  given  by  y  =  -h„h  —  ,  where  /r  (0)  S  1  and  the  length  scale  L  determines  how  fast 

the  undisturbed  fluid  depth^  is  changing.  Because  the  fluid  is  regarded  as  inviscid  we  may 
assume  a  velocity  potential  <f>  =  <f>(x,  y,  i )  exists.  In  terms  of  this  velocity  potential  the  equa¬ 
tion  of  motion  and  boundary  conditions  may  be  written: 

+  4>yy  “  0  for  <  y  <  T>(*.  7) 


Here  the  free  surface  is  given  by  y  -  rj(x,t),  where  /,  of  course,  is  (he  time.  The  first  equation 
represents  mass  conservation  at  each  point  in  the  fluid,  the  second  equation  expresses  the  fact 

•Such  questions,  while  of  general  interest,  arc  particularly  important  in  the  study  of  the  singularity  as  h  —  0 
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(hat  the  bottom,  y 


-Ml  4 


,  is  impenetrable;  here  h'(£)  =  — .  The  last  two  equations  are 


the  boundary  conditions  at  the  free  surface  y  «  rj.  Here  g  is  the  acceleration  of  gravity.  The 
third  is  the  kinematic  condition  at  the  free  surface  and  the  fourth  is  the  Bernoulli  equation 
evaluated  at  the  free  surface.  Implicit  in  the  derivation  of  the  latter  is  the  assumption  that  the 
pressure  is  constant  at  the  surface;  hence  we  are  ignoring  any  surface  tension  and  the  effect  of 
any  fluid  above  the  water. 


To  complete  our  description  of  the  fluid  motion  we  need  an  initial  or  boundary  condition 
which,  of  course,  produces  the  wave  motion  of  the  free  surface.  Fortunately,  we  do  not  need 
an  explicit  condition  in  order  to  derive  the  variable  coefficient  KdV  equation.  It  will  be 
sufficient  to  assume  that  this  condition  introduces  an  amplitude  a  and  characteristic  length  or 


wavelength  X  into  the  problem;  for  example,  a  typical  initial  condition  might 


icieristic  length  or 
be  ri  •«  a  at 


i  =  0.  The  initial  or  boundary  condition  will  only  enter  the  problem  through  a  and  X;*  these 
will  be  used  for  scaling  purposes  later. 


We  now  nondimensionalize  the  equations  of  motion  as  follows 

x  -  Xor 

y  “  h0y 

'"vfe' 

V  - 

4>  -  ek^/gFo  <t> 


where  x,  y,  /,  t).  <f>  are  the  new  nondimensional  variables  and  e  =  The  above  equations 


become 


h„ 


l*24>XX  +  4>yy  -  0  ~h(hx)  OOl| 

-  —n2&4>xh'(tix)  on  y  —  —  A(xft) 


4>y  “  +  ««MJ 


-t 7«*,  + j**,  ■  2 


on  y  -  «tj, 


(1) 


where  m  =  — ,  6  = 

A  L 


We  shall  now  study  waves  whose  amplitudes  are  small,  i.e.,  e  «  I.  We  shall  also 
assume  that  the  water  is  shallow,  that  is,  that  the  lengths  of  the  waves  are  large  compared  to 
the  depth,  i.e.,  fi  «  1.  In  addition,  we  will  regard  the  bottom  as  slowly  varying,  that  is,  the 
lengths  of  the  waves  are  taken  to  be  short  compared  to  L,  i.e.,  ft  «  1.  In  terms  of  the  physi¬ 
cal  parameters  this  requires  a  «  h0  «  X  «  L.  At  this  stage  we  need  not  assume  any  par¬ 
ticular  ordering  of  c,  n,  and  ft;  in  this  we  differ  from  Johnson  (1973)  who  at  least  at  the  outset 
takes  t  -  0(j*2)  -  0(ft). 

*lt  should  be  mentioned  that  we  are  also  assuminj  that  the  initial  or  boundary  condition  is  such  that  to  lowest  order 
the  wavc(s)  move  to  the  right,  i.e.,  in  the  positive  x  direction 
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Because  we  are  taking  «,  /x2  and  8  small,  a  natural  approach  to  this  problem  is  to  seek 
perturbation  solutions 

V  -  17 (a,  f;  t,  8,  /x)  -  r?u  +  n2i)t  +  8 172  +  +  M28t)s 

+  82t?,,  +  /aJer}7  +  8«tjk  +  e2T?9 

+ . 

and 

<t>  —  y,  /;  c ,  8,  /x)  —  4-  /xV(  +  8^2  +  *^3  +  /xV 4  +  /x28$5 

+  82(^6  /*2*^> 7  ■*"  8c +  *V9 

+  ... , 

where  the  dots  indicate  terms  of  third  order  in  c,  8  and  /x2  and  t?, -t?,(x,  /)  only  and 
<t>,  —  <Mx,  y,  /)  only,  /  -  0,  1,  2,  3 .  Although  we  are  regarding  c,  8  and  /x2  as  arbi¬ 

trarily  small  these  expansions  are  seen  to  break  down  at  large  values  of  xand  t.  In  particular,  it 

can  be  shown  that  t?,  contains  terms  which  are  proportional  to  x;  thus  when  x  —  O  the 

.  H 

M  i?  term  is  just  as  large  as  the  t?0  term  and  we  no  longer  have  a  valid  asymptotic  expansion. 
Nonuniformities  such  as  these  are  well  known  and  techniques  have  been  devised  which  can 
render  the  above  expansion  valid  for  large  values  of  x  and  /.  One  of  these  techniques  is  called 
the  method  of  multiple  scales  or  the  mullitiming  technique;  see  Cole  (1968),  Van  Dyke 
(1974),  Nayfeh  (1974)  or  Leibovich  and  Seebass  (1974).  In  fact,  both  Cole  and  Leibovich  and 
Seebass  have  used  this  technique  to  derive  the  flat  bottom  or  constant  coefficient  version  of  the 
KdV  equation,  suggesting  that  this  is  a  natural  way  to  derive  the  variable  coefficient  equation. 
The  basic  approach  of  the  mullitiming  technique  is  to  regard  7?  -  t?(x,  r;  e,  8,  /x2)  as  a  func¬ 
tion  of  new  independent  variables 

r  =  ^  *(*2)  "  '■ 

X0  X , 

x,  =  /x2x, 

X;  =  8x,  (2) 

and 


*3  =  «x, 

i.e.,  7?  -  t?(t,  x0,  x | ,  x2.  x3;/x2,  8,  e).  Here  r  is  the  wave  coordinate  previously  introduced. 
This  allows  for  the  change  in  the  wave  speed  due  to  the  changing  depth;  in  theories  where 
there  is  no  change  in  the  propagation  speed  r  would  simply  be  r  -  x  -  1  or  /,'(x2)  =  x?  =  8x. 
In  the  multitiming  technique  r,  x„,  xt,  x2,  and  x3  are  all  regarded  as  independent  variables. 

When  x„  -  x  -  ol-j-l,  r  —0(1)  and  x2  —  0(1)  and  when  x„  -  x  -  o[-^r|,  x,  -  0(1)  and 


when  x0  -  x  -  Oj— J,  ”0(1).  With  this  assumption  we  then  go  through  the  perturbation 

procedure  just  as  before.  However,  we  require  that  our  expansion  be  valid  for  x  -  O  -V,  — . 
»  M2  « 


and  —  j;  it  turns  out  that  this  requirement  gives  the  dependence  of  t>„  on  xt,  x2,  x3  and  ulti¬ 
mately  the  variable  coefficient  KdV  equation. 
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One  molivation  for  this  technique  is  found  by  examining  the  failure  of  the  straight¬ 
forward  expansion.  That  expansion  is  seen  to  have  the  following  form 

1 7  “  n„  +  M2!-*  '  const.  +  /(*))  +  . . . 

where  both  i)0  and  /  are  bounded  in  x.  This  would  seem  to  suggest  that  the  exact  solution  for 
T/  can  be  written 

1)  =  7)(x,  /;  € ,  /a2,  8)  =  i)(x,  t,  fl2 x,  ix2t\  €,  /x2,  8). 

Similar  remarks  apply  for  €  and  8.  Another  more  intuitive  explanation  is  that  we  see  two 
simultaneous  changes  in  x.  The  first  is  due  to  the  fact  that  the  wave  carries  changes  in  the 
water  height  and  in  the  fluid  velocity.  By  definition  of  X  these  changes  occur  over  a  distance 
x  =  0(1).  In  a  coordinate  system  traveling  with  the  wave  these  are  seen  as  variations  with 
respect  to  t.  The  second  is  that  the  basic  wave  shape  is  being  slowly  modified  due  to  non¬ 
linearity,  dispersion  and  the  bottom  variation.  Thus,  it  would  seem  reasonable  to  express  t)  in 
the  above  manner.  Of  course,  mullitiming  techniques  are  well  established  in  mechanics  and 
applied  mathematics  and  these  comments  are  not  meant  to  be  a  complete  justification  of  the 
technique,  rather  they  are  meant  to  remind  the  reader  of  some  of  the  intuitive  ideas  connected 
with  this  method.  For  a  more  complete  discussion  of  multitiming  techniques  the  reader  is 
referred  to  the  above  mentioned  references. 


3.  Derivation  of  the  Variable  Coefficient  KdV  Equation 

In  our  derivation  we  first  need  to  expand  equations  (I)  in  powers  of  e,  8  and  n2.  In  the 
following  we  may  regard  e,  8  and  fi2  as  independent  and  of  any  relative  size.  However,  as  in 
any  multitiming  scheme  involving  more  that  one  small  parameter  there  will  be  difficulties  with 
this  viewpoint,  see,  for  example.  Chapter  IV  of  Leibovich  and  Seebass  (1974).  Because  these 
are  most  easily  seen  by  an  analysis  of  higher  order  terms  not  discussed  below,  we  shall  delay 
the  discussion  of  this  until  Section  4.  It  will  be  noted  at  this  stage  that  these  difficulties  are  to 
be  expected  in  any  multitiming  approach  and  of  course  are  implicit  in  Johnson’s  (1973)  work. 


We  shall  first  apply  the  basic  assumption  of  mullitiming,  namely  that  the  independent 
variables  are  given  by  (2).  We  do  this  in  a  stepwise  fashion;  first  we  change  variable  from  x 

and  t  to  x  and  t,  where  t  =—  j?(8x)  -  t  and  x  =  x.  The  x  and  r  derivatives  can  then  be 
transformed  using  the  chain  rule 


and 


JL 

dt 

JL 

dx 


dx 2 


JL 

dr 


JL 

dx 


dxdr 


,  92  ,  j. ■■  9 


_  jjs 
dt 

and  r;  again  by  the  chain  rule 


where  #(f )  =  — .  We  further  transform  x  and  r  to  the  independent  variables  x()l  Xj 


9 1 

_9_ 

dx 


JL 

9r 


9 


9 


9 


*£  +  ^  +  ‘*,'£7  +  8977  +  <-k1 


5 


MARK  S.  CRAMER 


and 


dx2 


-  +  2g—‘—  +  — 

dr2  dx$T  8x02 


32 


+  5 


*  * +  u  *  2j£t, 


+  2 m2  i 


a2  .  a2 


3X|3r  dX|dx0 


+  2,1*^  +  »’ 


dx3dr  Sxjdx0 
+  0(#i4,  /x2S,  /i2«,  82,  Se,  t2), 

here,  of  course,  g  =  g(x2)  and  in  (1)  we  will  take  h  —  /»(x&)  —  h(x 2).  Equations  (1)  there¬ 
fore  become 

<t>yy  =  ~^U2<t>rr  +  lg4>x„r  +  4>x„x)  +  2 g<t>X}T  +  W>x„x^ 

-  2nAlg  <*>X|T  +  4>X  X)]  -  +  ♦xAI 

+  Olp.6,  ja48,  f44€,  p.282,  fi.28€,  fi2«2l,  (2a) 

<t>y  -  -/u28(£^T  +  <f>x  )  h'(x 2)  +  o(m28)  on  y  *  ~h(x2)  (2b) 

=  — m*»It  +  €/*2(gqT  +  ■»)*  )  (g$r  +  4>x  )  +  oitfjL2)  on  y  «  er),  (2c) 


and 


and 


and 


-TJ  «  ~<£r  +  —  €  (g<ftT  +  0X))2 


+  *f*26x,(i4>T  +  4>x„)  +  +  4>x,) 


+  €2<f>X}(g<l>T  +  lf>xJ  +  0(«/*4,  82€,  €3) 

+  4-^7  4>2,  on y  -  €T).  (2d) 

l  H 

At  the  next  stage  of  the  analysis  we  would  need  to  expand  <t>  in  the  two  boundary  conditions  at 
y  -  «tj,  eqns  (2)c,d  on  y  -  erj  we  may  write 

<f>  -  <£(x,  t.  y)  *  <t>(x,  l,  er)). 

Hence  for  small «  this  may  be  expanded  in  a  Taylor  series  about  y  —  0: 

<t>  -  <b  (x,  /,  0)  +  «tj  +y  (x,  t,  0)  +  . . . 
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In  like  manner  the  derivatives  *T.  *v  4>X],  *Xj,  4>X},  4>y  appearing  in  (2)c,d  may  also  be 
expanded 

*t(«tj)  -  *T  (0)  4-  er)  <f>yr(0)  +  ...  , 

+  er) 

*x,  («Tj)  -  *X|  (0)  +  et)  4>yXi  (0)  +  . .  .  , 

<t>x2(*v)  -  4x}(0)  +  er)  <l>yX2( 0)  +  ... ,  (3) 

*x3 («>)  -  +  er)  **x3(0)  +  •  •  ■ . 

and 

*y(«T })  “  *,(0)  +  er)  <f>yy( 0)  4-  . . .  , 

Ordinarily  we  would  simply  substitute  these  expansions  in  the  boundary  conditions  (2)c,  d. 
However,  to  save  space  we  shall  anticipate  a  result  which  comes  out  or  the  analysis  presented 
later.  This  is  that  4>y(r ,,  x0,  X\,  ....  y)  —  0(ji2)  at  most.  We  now  use  this  result  and  substi¬ 
tute  expansions  (3)  in  equations  (2)c  and  (2)d;  these  boundary  conditions  then  become 

*v  *  ~er)  -  n2t)T  +  efi2(gr)r  4-  r)x  )  ( g<f>r  +  <f>x)  4-  oiep2),  on  y  «  0,  (4) 

and 

-r)  -  -*T  -er )<f>yT  4-  y«(£*T  4-  4>x)2  +  e/x2*Xl(£*T  +  <f>x) 

+  eS  *,j(£*t  4-  <f>x  )  4-  «Vx3(i^r  +  *x„) 

4-  0(e2fi2,  efi4,  S2e,  e})  4-  ~  *2, 

on  y  -  0.  (5) 


Because  we  are  anticipating  the  result  < py  -  0(jt2)  we  see  that  the  main  effect  due  to  the  expan¬ 
sions  (3)  is  the  addition  of  a  -er) <f>yy  term  to  (2c)  and  -er)4>yT  term  to  (2d)  and  the  transfer  of 
the  boundary  condition  from  the  moving  surface  y  —  er)(x,  t )  to  the  stationary  surface  y  —  0, 
at  least  to  the  order  considered  here.  The  surface  elevation  t)  and  the  potential  *  are  now 
given  by  the  equations  (2a),  (2b),  (4)  and  (5). 

Finally  we  need  to  expand 

■t)  -=  t)(r,  x0,  x 1,  x2,  x3\  e,  fi2,  8), 

and 


*  -  *(r,  X0,  Xj,  x2,  X3;  €,  /x2,  8) 

in  expansions  of  powers  of  e,  n2  and  8,  i.e., 

V  "  Vo  +  +  fri?2  4-  er)}  4-  /x4 r)A  4-  /x28t)5  4-  8*»}t  4-  fi2er)7  4-  8er)8  4-  «*!)9  4-  . . .  . 

^■^o+  n2<t>\  4-  8*2  4"  «*  j  4-  /x4*4  +  M2®*5  +  82*6  4-  /x2«*7  4-  8e*g  4-  c2*?  4-  ...  . 

Here  the  dots  refer  to  terms  of  third  order  in  /a2,  8  and  c.  When  these  expansions  are  substi- 
t  ited  in  equations  (2a),  (2b),  (4)  and  (5),  and  terms  of  like  order  in  #i2,  8,  «,  /x4,  fi2b,  ... 
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are  equated,  we  obtain  a  system  of  equations  which  need  to  be  solved;  these  are  listed  in  the 
Appendix.  We  now  proceed  to  solve  these  keeping  in  mind  that  we  want  to  eliminate  any  secu¬ 
lar,  i.e.,  proportional  to  x0,  terms  which  may  occur  in  tj,  or  <f>r 


From  equations  (Al)  we  see  that  d>0  =  B0(t,  x0,  x\,  x2,  x3)  only.  We  will  also  require 
that  the  wave(s)  move  in  the  positive  x  direction  only  and  will  therefore  assume  that  B0  is 
independent  of  x0,  i.e., 

<t>0  =  B0(r,  X,,  x2,  x3), 
and 


Vo  Bot. 

We  now  consider  equations  (A2);  these  now  read 

$  \yy  *  8  B011t 

<t>\y  =  0  on  y  =  —h, 
and 


<t>\y  “  ~Bor T  on  y  =  0. 

The  first  equation  implies 

4>\y  -  -i2  B0„  y  +  A{(t,  x0,  xu  x2,  x}) 
Substitution  in  the  boundary  condition  at  y  =  —  h  yields  A\: 

A  i  =  ~g2h  B0„, 

hence 


<t>\y  =  -g2  B0TT(y  +  h)  in  general. 

We  also  need  to  satisfy  the  boundary  condition  at  y  =  0;  this  then  requires  that 

g2  =  _J_ 
h(x2) 


or  upon  integration 


g(x2)  “  JQ  2  h  v’  (p)  dp. 


where  we  have  chosen  the  positive  sign  for  g.  This  of  course  agrees  with  the  result  presented 
earlier  for  t.  Thus,  our  results  for  d>\  and  tji  are 


Bott  +  B |(t,  X0,  X!,  Xj,  X}), 


(6) 


and 


Vi  “  B\t 


We  now  consider  equations  (A3).  Integration  and  application  of  the  boundary  conditions 
at  y  —  -h  and  y  -  0  show  that  <f>2v  =  0  and  therefore 

<t> 2  -  B 2(t,  x„,  x,,  x2,  x3), 


and  from  the  last  equation  in  (A3) 
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We  now  consider  equalions  (A4).  As  in  equations  (A3)  we  can  easily  show  that 
<t> 3  =  B)(t,  x0,  x,.  x2,  *3) 


T)3  -  BiT  -  —  7)1 


When  we  use  the  fact  that  <j>ox  =  0  and  our  results  for  <t>\,  equations  (A5)  become 

2 

$4 yy  ~  ~  B]T7  +  ~j=  BixaT  4"  +  ^2  2  ^  Vottt  —  Vox  ,• 


04V.  =  0  on  y  =  —h 
<tuy  =  -fl|TT  on  >  =  0. 

If  we  integrate  the  first  equation  we  find 

I' 

ir +  ^2~  t)oitt 


Vox |  T  ■*<>»  •*! »  -*2»  ■^3^* 


The  boundary  condition  ai  y  =  —h  determines  /f4  as 

A 4  =  -2Va  T}ttV|  -  Y  7)tlTrT  -  h  —  BUt  +  -i-  Bix  +  Bix  x 


The  boundary  condition  at  y  =  0  gives  us  the  equation  for  B\\ 

A  4  =  ~B\jt 


B\xnT  "I"  B\Vi  —  2x/h  7)nXi  +  J  VlITTT 

We  note  that  a  particular  solution  for  B\  is 


Bi  =  -x„^f-  /  Ktdr 


t)i  —  B\,  —  -xa~-  AT | , 

where  AC,  =  K\(r,  xt,  x2,  x})  =  2-Jh  tj„m  +  -j  T)„rrr.  Thus,  t?|  grows  linearly  with  x„  and 
when  x„  —  x  —  o[  — ^-1  our  expansion  scheme  breaks  down,  i.e.,  /xJrj|  -  0(t)„),  unless  Af,  =  0, 


9 


MARK  S.  CRAMER 


,  yfh  n 

VoXf  ^  j  VoTTT  “ 


(7) 


1 


Thus,  in  order  to  avoid  the  breakdown  of  our  expansion  when  x  -  0|— yj  we  shall  assume  that 

the  x |  behavior  of  t)b  is  given  by  (7).  We  note  that  there  is  no  assumption  here  on  B\,  except, 
of  course,  that  BiT  -  r),  be  such  that  our  expansion  be  valid  for  large  values  of  x.  A  more 
careful  analysis  of  the  homogeneous  solution  for  B\  shows  that  both  right  and  left  moving 
waves  are  possible.  Had  we  assumed  Bix<>  =  0  or  r)tXo  =  0  as  does  Johnson  (1973),  this  would 
not  be  the  case. 


We  now  consider  equations  (A6)  which  now  read 


<t>5yy  ”  - 


h  Bl7T  +  yfh  Blx“r  +  B2x“x" 
h' 


+  y  h  2  h'  Vo  -  t)0 


<t>sy  =  Vo  on  y  -  -h 

<t>$y  =  -&2TT  on  y  -  0. 
Integration  of  the  first  equation  yields 

<t>iy  “  ~ 


i  B2rr  +  -Jp  B2Xor  +  B^ 


y  +  ~  h  2  h'i)0y 


7 1 


Vox.y  +  ^5(r,  X0.  X | ,  x2,  x}) 


The  first  boundary  condition  yields  A 5 

1  A' 

^5 


1  JL 

2  VA 


7}0  -  2yfh  Vox 2  -  \  B2tt  +  B2x„r  +  B2Xiixo  h. 


The  second  boundary  condition  yields  our  equation  for  B2. 

A  5  “  -B2l1 


7  D  i  D  — _ L 

JJj  a2v„  fj 


27^  Vox2  +  2  xfh  Vo 


A  particular  solution  for  B2  is 


or 


B2=~x 0^~  f  K2  dr 


V2  "  B2r  -  -X0~-  Kl • 


where  K2  =  -j- 
h 


«  rr  ,1  h' 

2  V/r  t,0X2  +  J^Vo 


.  Thus,  when  -  x  -  0|  — 

8 


,  8 172  "  0(»?o)  unless  Aj 


=  0.  Thus,  in  order  to  avoid  a  breakdown  in  our  expansion  when  x  »  0  —  we  shall  require 
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If  we  consider  equations  (A7)  we  can  easily  show  that  ^6  “  B4(t,  x0,  jt|.  x2,  *j)  and 
B(n- 

When  our  results  for  <t>\  and  4>}  are  substituted  in  (A8)  we  have 


12  2 
tlyy  ~  ~  Byr  +  Bjv  +  Bix0x0  -  Voxy 


fay  —  0  on  y  —  —h 
fay  -  VoVot  ~  Birr  on  y  -  0. 
Integration  of  the  first  yields 


*ly  “  ~[i  +  +  y~7hVox>y 

+  Aj(t,  Xg,  X\.  X2,  x3). 

Substitution  into  the  first  boundary  conditions  shows  that  A4  is  given  by 


A-l  -  -  ~  Birr  +  B3xBr  +  B ^  h  -  2 -Jh  T)„3 


The  second  boundary  condition  yields  our  equation  for  By 

3 

Aj-  -Byr  +  ~fj  VoVor 


~T T  Bix„r  +  B)Xoxo  **  ^  2yfh  TJaxj  +  h  VoVor  ■ 


A  particular  solution  for  B}  is 


B,--x0^f  K}  dr 


By  -  Ky 


where  tf3  =  ^2-fh  n«*5  +  |  t>0tjot|.  Thus,  ijj  -  By  -  \  ~  vi  grows  linearly  with  *0  and 

when  x„  -  x  -  o|-j-J  our  expansion  breaks  down,  i.e.,  tr\i  -  0(rj„),  unless  K}  =  0.  Thus  we 
shall  require  that  A3  =  0,  that  is, 


3  ~~ 

Vox }+  J  h  1  Wot  =  0. 


If  we  continue  on  to  equations  (A9)  and  (A10)  we  see  that 
fa-  Bt(r,  xot  x]t  x2.  x3)  only, 

fa  -  B4(t,  x0,  X\,  x2.  x3)  only. 
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The  functions  v«  and  179  are  of  course  given  by  the  last  equations  in  (A9)  and  (A10). 


This  completes  our  analysis  of  the  perturbation  equations.  In  order  that  our  expansion  be 


valid  up  to  x  -  Ol  -~r 


or  O 


— |  or  O 
€ 


we  see  that  we  need  to  modify  the  behavior  of  v0  at 


those  large  distances;  this  is  the  expected  result  of  any  multitiming  scheme.  The  necessary 
modification  to  rj0  are  given  by  equations  (7),  (8),  and  (9)  and  are  rewritten  here  for  conveni¬ 
ence 


and 


dVo  =  d3Vo 

dXj  6 

dy0  _ j_  h'(x2) 

dx2  4  h(x2)  1*° 


dy0 

dx3 


3 

2 


dy0 

dr 


We  can  show  that  t)0  must  satisfy  the  variable  coefficient  KdV  equation  by  "reconstituting”  the 
above  equations,  that  is,  by  writing 


dx 


dy0 

dX0 


+  fl‘ 


jh)o 

dx. 


+  82^ 

(r.x0.xj.x,  OX2 


+  e  aT>0 

(r.jr„,8-|,xj  dx3  jr,x0.x,.x2 


This,  of  course,  just  follows  from  the  chain  rule  and  has  already  been  used  to  transform  from 
x,  r  coordinates  to  the  r,  x0,  X|,  x2,  x3  coordinates.  Thus, 


dy0 

dx 


2yfh 

M  ^  VlOTTT 


5  W 
4  h Vo 


3 


£  2  * 


Vo  Vo 


or 

Vox  VoTTT  +  1  ^oVoT  “  0>  (10) 

where  we  would  now  regard  h  =  hi fix).  Once  we  allow  for  differences  in  notation  we  see  that 
Equation  (10)  is  exactly  equation  (14)  of  Johnson  (1973)  or  equation  (l)  of  Ostrovsky  and 
Pelinovskiy  (1970). 


4.  Discussion 

We  have  derived  the  variable  coefficient  KdV  equation  (10)  by  standard  multitiming  pro¬ 
cedures.  This  gives  the  correct  behavior  of  Vo  over  distances  and  times  of  order  —  or  ~  or 

«  1 


whichever  gives  the  largest  distance.  The  behavior  for  larger  distances,  say  x  -  o|-y 

aI  1  I  .. _ 1  j _ . _ _ _ 1 _ •  r  1  ‘  ■ _ _ 1 _ .  _  ■  _  ■  ■  l  .1  . 


or 


Ol — r  ,  would  require  an  analysis  of  higher  order  terms  and  would  show  that  we  need  to 

W  ) 

further  correct  t )„,  Thus,  some  care  should  be  exercised  when  one  of  the  parameters  is  very 
much  less  than  the  other  two.  For  example,  if  8  «  e  or  fi 2.  say  S  **  cJ,  then  (10)  is  not 

expected  to  give  an  accurate  description  of  the  flow  for  | 


include  the  corrections  necessary  at  x  —  0 


We  would  then  need  to 


-yl.  In  this  case  (10)  would  reduce  to 
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u2  3 

Vox  r  VoTTT  +  T  e  VoV  or  “ 


and  would  onlv  be  valid  for  x  -  0  — ,  — r|.  One  obvious  way  to 

I*  P  \ 

t.  8  and  m :  are  all  much  larger  than  any  quadratic  term  involving^, 
ensure  that  the  only  corrections  needed  for  x  -  of— ,  -j-j  are 


M 


8 


(10).  Other  than  this  the  relative  sizes  of  e,  pr  and  8  are  arbitrary. 


avoid  this  is  to  require  that 
e  and  p2.  This  would  then 
those  contained  in  equation 


Another  point  of  interest  is  that  the  analysis  presented  here  does  not  preclude  the  possi¬ 
bility  of  reflected  waves  in  tj(.  t)2,  and  173.  A  careful  analysis  of  these  terms  shows  that  both 
right  and  left  moving  waves  are  possible  in  these  solutions  although  the  analysis  presented  here 
is  only  expected  to  give  us  the  behavior  of  tj0  and  not  that  of  the  higher  order  terms. 

In  conclusion,  the  analysis  presented  here  is  somewhat  more  general  than  that  of  Johnson 
(1973).  The  derivation  presented  here  is  in  terms  of  a  standard  technique  of  singular  perturba¬ 
tion  theory,  namely  multitiming.  The  variable  coefficient  KdV  equation  (10)  gives  an  accurate 
description  of  the  flow  for  distances  of  order  of  the  largest  of  c-1,  p~2  and  8~'.  Except  for  the 
difficulties  mentioned  earlier  in  this  section  there  are  no  other  restrictions  of  the  size  of  e,  p2 
and  8  relative  to  each  other. 


This  work  was  performed  at  the  Naval  Research  Laboratory  while  the  author  participated 
in  the  ASEE  Summer  Faculty  Research  Program.  Many  useful  discussions  with  James  Witting 
formed  part  of  the  basis  of  this  work;  his  support  is  gratefully  acknowledged. 
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Appendix  A 

Equations  Resulting  from  Perturbation  Expansions 


The  perturbation  procedure  outlined  in  Section  3  leads  to  the  following  system  of  equa- 


ions: 

0 oyy  “  0 
0ov  “  0 

0„>“O  j 
Vo  “  0ot| 

01 »  “  ~^20orr  +  2 g4>ox„r  +  0<w„v„l 

0I.V  »  0 
0ly  ”  *>otJ 
V\  * “  0  It  I 

02v.  “  0 

02y  “  0 

02V-O  j 

V2~  02rJ 

03yy  ”  0 
0JV-O 

03>  “  0 

1fl3-03r  “  J  (0or *  +00*,,)  2 

04w  “  _li20trr  +  2j0Ur  +  0|v„v„  +  2#0, 

04v  -  0 
04y  “  ~  V  !t| 

V4  “  04r  I 


on  y  —  -h  (All 

on  y  -  0 

on  y  —  —h  ( A2) 

on  y  “  0 

on  y  -  —h  <A3) 

on  y  —  0 

on  y  «  -h  <A4) 

on  y  —  0 

,r  +  20ov„v,) 

on  y  -  (A5) 

on  y  -  0 


05W  ”  “li202rr  +  2g$lx,  +  <t>2x„x„  +  <T0or  +  2g0OAjT  +  20,,,,,) 

05v  “  “(i0or  +  0<w„l  A'(*2>  on  y  -  -/t  ( A6) 
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*5y  “  _t)2t 
VS  “  *5, 
*6„  "  0 
*6,  “  0 

*6,-0  | 
n6  -  *6,  J 


on  y  «■  0 

on  y  —  -h 

(A7 

on  >  —  0 


*7„  -  +  2j*J,ir  +  *3^.,,  +  lg+ax,,  4-  2*nj(j,  ] 

*7,-0 

*7,  —  n#i*l„  ^3,  ^iViiT  9#,,^)  *«jt#) 

-  7)7  -  *7r  -  T)o*l,r  +  *or*lr  +  *«„*U„  +  "j"  *D„*or  +  "j  *It*ox„ 

+  *<u,**or  +  *«,*«„  +  -J  *ly 


on  y 


on  y 


-h 

(A8) 

0 


** yy  “  0 

*»>  —  0 

*«,  -  0 
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In  these  equations  we  have  anticipated  and  used  the  result  that  *„  *  *0(y)  and  in  equations 
(A9)  and  (AIO)  we  have  used  the  result  that  *2,  -  *3,  =  0. 


